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Abstract 



We introduce a class of probability measures whose densities near infinity are mixtures 

of Pareto distributions. This class can be characterized by the Fourier transform which has 

a power series expansion including real powers, not only integer powers. This class includes 

Q> ■ stable distributions in probability and also non-commutative probability theories. We 

P^ . also characterize the class in terms of the Cauchy-Stieltjes transform and the Voiculescu 

fH I transform. If the stability index is greater than one, stable distributions in probability 

"t^ ' theory do not belong to that class, while they do in non-commutative probability. 
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en ■ 1 Introduction 

O I Probability measures with power-law behaviour have found their applications in a variety of 
phenomena such as the energy spectrum of fluid [6l [TT] , the distribution of dark matters [18] 
and deformed Gaussian distributions in Tsallis statistics [23] • Readers interested in further 
. information on power laws can find reviews such as [T^]. A basic power law in applications 
i^ ! is the Pareto distribution whose density is expressed as cx~"~^ {x > R > 0) ioi a > 0. In 
p^ \ particular, a noise is called a 1// noise, or more precisely 1//"+^ noise, if the frequency / follows 
a Pareto distribution with parameter a. The case — 1 < a < also appears in applications. In 
this case, the integral J^ f^^^^df diverges. In practice, this does not cause a problem since 
we usually focus on a finite interval [i?, R'] . However, we only consider a > in this paper for 
mathematical simplicity. There are a large number of researches on the 1// noise. For instance, 
an explanation of the origin of 1// noise can be found in [3]. 

In this paper, we introduce a class of probability measures which are mixtures of Pareto 
distributions. More precisely, let S be an additive sub-semigroup of [0, oo) such that: 

(SI) nas- 

(S2) there is a constant c > such that ^{S r\[n,n + 1)) < c"+^ for any n G N. 
Then we introduce the following class of probability measures. 



Definition 1.1. A probability measure /i is said to be in Vs if the following condition is 
satisfied: there exist {ap}p^s,i3>o C C and r, i? > such that 

iJ'\\x\>R{dx) = ^ Im I a/3 ( - j ] dx, (1.1) 

l«/3|<r^ (1.2) 

r < c~^R, (1.3) 

where c is the constant in the condition (S2). The function {-)'^'^^ is understood to be 
gi(/3+i)7r|^|-/3-i £qj, ^ ^ _^_ Under the conditions (11.21) and (II. 3p . the series in (II. ip is ab- 
solutely convergent (see Proposition 12. 4p . 

Let J-'^{z) = f^e^^^fi{dx) and G^{z) = J^j^^{dx) be the Fourier transform and the 
Cauchy-Stieltjes transform (or Stieltjes transform for short) of a probability measure n, re- 
spectively. The main results of this paper are as follows (see Theorems 12.61 and 13. 3|) . The 
following are equivalent: 

(1) /^ e Vs; 

(2) There exist (c^)^^^ C C and A > such that \c^\ < A^ for any 7 G 5, Cq = 1 and 

(3) There exist {d^)^(zs C C and A > such that \d^\ < A'^' for any ■y & S, do = 1 and 

^^(^) = E;7iT' ^eC_, 1^1 >cA (1.5) 

-yes 

The above equivalence generalizes a result of p] which corresponds to the case S* = N. The 
power z'^ is defined as the principal value in C\(— oo,0]. The coefficients c^ and d^ in fact 
coincide. In Subsection 13.21 we moreover characterize Vs in terms of the reciprocal Cauchy 
transform and the Voiculescu transform. We consider G^ on C_ := {2 G C : Imz < 0} rather 
than on the upper half-plane. This is because the lower half-plane for G^ plays the role that 
the positive real line (0, 00) does for J-"^. This correspondence is clarified more in Section [31 

If we calculate the Fourier transform of the Pareto distribution restricted to [0, 00), a loga- 
rithmic term z^ logz or other singular terms may appear (see Proposition 12. 5p . To cancel such 
a term, we consider Pareto distributions on the negative real line, not only on the positive real 
line. However, sometimes we do not have to consider the negative part. In Subsection 12.31 we 
clarify when Pareto distributions on the negative real line can be removed, using a Diophantine 
approximation. 

The idea for the equivalence between (1) and (2) comes from stable distributions [2S]- The 
distribution /x of a random variable X is said to be stable if for any a,b > 0, there are c > 
and d ^M. such that aX + bX' has the same distribution as cX + d, where X' is an independent 
copy of X. If c? is for any a, b, the distribution n is said to be strictly stable. Except for the 
Gaussian case, the density of a stable distribution shows a power-law tail cx~°'~^ as x — )■ 00, 
a G (0,2). The parameter a is called an index of stability. A stable distribution has two 
important aspects: the distribution of a self-similar Levy process [2D1; the limit distribution of 
the sum of i.i.d. random variables |[7j. 



A stable distribution /Xq, with stability index a has a Fourier transform of the form 

^ r^N /^». ., x_ /exp(27^-c(l-z/3tan(f)-sign(z))|z|"), a 7^ 1, 
Jm [exp (27z-c(1 + «/^- -(log|2|)sign(2;))|2;|j , a = 1, 

where c > 0, 7 G M and /3 G [—1, 1] [2D]- We do not consider the case a = 1 and /3 7^ 0. Then 
we can use a different parametrization: for z > 0, 

J-^Jz)=exp(z7;2 + ^'^6z'^), (1.7) 

where h G C\{0} satisfies arg6 G [(1 — a)7r,7r] if < a < 1 and argfo G [0, (2 — a)7r] if 1 < a < 2. 
It suffices to consider 2; > since information on z < can be recovered from the complex 
conjugate. The Fourier transform J^f^^iz) can be expanded in the form 

00 .^ 

which is a special case of (11.41) with Sa '■= {fn + na : m,n E N}. However, the coefficients c^ 
are bounded as claimed in (2) if and only if a G (0, 1]. 

Since the coefficient c^ (= d^) in (ll.4p has an analogy with the moments of a probability 
measure, we call c^ a 7-complex moment. This is why we avoid a logarithmic term in the 
Fourier transform: the concept of 7-complex moment can be defined thanks to the absence of 
a logarithmic term. 

The idea for the equivalence between (1) and (3) comes from stable distributions in non- 
commutative probability, or more specifically, in free, monotone and Boolean probability the- 
ories [211 [161 EI]- The reader is referred to Section |3] for details and we just mention here 
a remarkable difference between probability and non- commutative probability theories: for 
a G (1,2], the estimate |c/3| < A'^ in (II. 4p does not hold in probability theory, while it does in 
free, monotone and Boolean probability theories. 

2 The characterization in terms of the Fourier transform 

2.1 Generalized power series expansions 

In this paper, N denotes the set of natural numbers {0, 1,2,3, ■ ■ ■}. z^^ := e^^°^^ for /3 G M 
denotes an analytic function in C\(— oo,0], where logz is defined so that — tt < Im(logz) < vr. 

Remark 2.1. (1) The following calculations are correct for z G C\[0, 00): 

(i) z^'z^ = z°+'' for a, /3 G M, 

(n) z"° = (z°)'" for a G M, ra G Z. 
However, (wz)"' 7^ w°'z°' and z°''^ 7^ (z")^ in general. 

(2) It is sometimes useful to denote the expansion (ll.Sp as 

T(z) = y -f ^I!1L^™+- ;, > 0. (2.1) 

However, the coefficients Cm^n may not be unique. This is because mi + n\a = 1712 + n2a can 
occur for distinct {1711,121) and (m2, 112). To avoid this non- uniqueness, we introduce summations 
over a set such as J^s^Sa rls+i) ^'^ ■ Then coefficients C/3 are unique. 



Definition 2.2. A series of the form 

7es 

is called a generalized power series. 

We show some examples of the set S satisfying the properties (SI) and (S2). 

Example 2.3. (1) Let a > and Sa '■= {m + na : m,n E N}. Then Sa satisfies (SI) 
and (S2). The condition (S2) can be checked as follows. If A; < m + na < k + 1, then 
k — m < na < k — m + 1. For each m E N, the number of possible n's is at most [-] + 1. 
Therefore, \Sa n [A;, A; + 1)| < ([i] + 1)A;. 

(2) More generally, let < ai < ■ ■ ■ < a^ and Sai,---,ap '■= {^o + ^i^i + ■■■ + n^Op : 
no, ■ ■ ■ ,np E N}. Then Sai,---,ap satisfies (SI) and (S2): it can be proved that a constant c > 
exists such that \Sai,---,ap D [A;, A; + 1)| < ck^. 

We prove a basic property on convergence of generalized power series. 

Proposition 2.4. Let {a^)^(zs be a sequence of non-negative real numbers. Then the generalized 
power series '^^(^^ ^i^"' converges in (0, e) for an e > if and only if there is an A > such 
thata^ < A^+^ 

Proof. We assume that the generalized power series Yl^es ^i^^ converges in (0, e) with e G 
(0, 1). We note that the estimate 



7G5 n=0 ■y£Sn[n,n+l) n=0 \ 7G5n[n,n+l) 



ay I z 



holds. This implies, from the usual theory of series, that aC > e ^ exists such that '^y^sn\n n+i) ^i — 
C"+^ for any n G N. In particular, a^ < C"+^ < C^+^ for 7 G 5 n [n, n + 1). 

Conversely, let us assume that there is an A > such that a^ < A'^^^. We can moreover 
assume that A > I. Then E^gs^T^^ ^ Er=o E^gsn^n+i) ^^^^^^ < cA^ J27=oic^^T is 
convergent in (0, {cA)~^), where c > is a constant such that IS* fl [n, n + 1)| < c"'~^^. D 

2.2 The main result on the Fourier transform 

Let [P] denote the integer part of /3, that is, [/3] is the largest integer which is not larger than 
/3. First we calculate the Fourier transform of a Pareto distribution. 

Proposition 2.5. (1) For P > 1, z > and R > 0, the following expansion holds. 

j'OO 

R^ / e'^^'x-^-^dx 
Jr 

= L^...(^ + 1^,)^ +^...(^-|^|+2)<^^> (2-2) 
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where Ci(/3) := J^ e^^x ^'^^^^ "^dx. fp{z) is defined as follows. ■ 






//[/S] = 1, i/ie /irsi summation of 112. 2\} is understood to be and /3 ■ ■ ■ (/3— [/3] + 2) is understood 
to be 1. If (3 ^ N, fp^z) can be included in the second summation of h2.2^) : however fp^z) is 
still exceptional for /3 ^ N since it may contain a logarithmic term under the limit /3 — )■ cxd. 
(2) For 0</3<l, 2>0 and R> 0, the following expansion holds. 



d^'x-^-Hx = C2{l3){Rzf + J2' \uB-k) ^ (2.3) 

-' fc=0 '^'^ ' 

where C2(/3) := J^ e'^x-^^-^dx + J2T=o wb)- 

Proof. (1) Let us assume that /3 > 2. For z > and i? > 0, we have 

/•OO /"OO 

R^ / e'^'^x-^-^dx = {Rzf / e'^'x-^-^dx. (2.4) 

J R Jrz 

With integration by parts, /^^ e*^x^^^^(ix can be computed more. 

/•OO 

Jrz 



/3 



P Jrz 



_ JRz) ,^^ i{Rz) + j^^ I 1"^ p'^r-l^+^dT (2.5) 

" /3 /3(/3-l) m-l)JRz 

^/3---(/3-fc + l) +/3---(/3-[/3] + 2)A/ "^ 
To expand the last integral in terms of powers of z, we decompose it into two parts: 
gix^-/3+[/3]-2^^ = / e'^x-^^^^^-^dx + / e^^a;-^+[^l-2rfa; 

= Ci(/3) + Yl V\ f ^'"'^^^^^"^dx + t f x-^+[^l-Mx .2 6) 

where gp{z) = i ''^%l'l^]^~^ if /3 ^ N and g(s{z) = -t\og{Rz) if /3 G N. (I22D follows from the 
above relations f l^ -f ETB]) . 



Let us consider the case 1 < /3 < 2, or equivalently, [/3] = 1. Using the same method, one 
can check that 

Therefore, the relation (12 ■2p is still true. 

(2) If < /3 < 1, the same method is applicable and then fl2.3p follows. D 

Therefore, a logarithmic term appears in the Fourier transform of x~^~^,x > i? if /3 G N. 
Even in the case /3 ^ N, a singularity can appear if we consider asymptotic behavior as /3 — )• oo. 
We remove this singularity by taking an appropriate function supported on the negative real 
line. Then we can characterize the Fourier transforms which can be expanded by z'^ for 7 G S. 

Theorem 2.6. Let ft be a probability measure. Then fi E Vs if and only if there exist c^ G C 
and y4 > such that: 

(1) \c^\ < A^ for any 7 G S", 7 > and cq = 1; 

(2) 7,{z) = E,e5 rETT)^'^' for z > 0. 

Remark 2.7. To generalize the expansion J-'fj.{z) = X^^o ^^^^^! (^^)"' fo^ ^ compactly supported 
H, we use the factor p, -*, -^^ for /i G Vs- In Subsection 13. 2^ this generalization turns out to be 
relevant. 

Proof. We prepare a notation for global behaviour of functions. For complex-valued functions 
f,g on S X N X [0, 00), f -< g means that there is a constant C > 0, independent of {f3,n, z), 
such that \f{/3,n, z)\ < C\g((3,n, z)\ for any (/3,n, z). 

Let us assume that /i G Vs- By definition, fi\\x\>R can be written as 



fJ'llxlyRi.dx) = y^ Im a/3 - ] dx, (2.7) 




where |a^| < r^ for an r G (0,c ^R). Let us decompose /i into three parts: /i = /i_ + /xq + 
/i+, where /i_ := /i|i,.<_i{, /io := /i||x|<_R and /i+ := fi\x>R- Since /iq is compactly supported, 
the Fourier transform of yUo can be expanded in a series Yl'^=o '""^"■^ z"^:" with the estimate 
|m„(/io)|<i?". 

(Step 1) We start from the cancellation of singular terms of J-"^^ and J^fj,_. Let us calculate 
the contribution of the Fourier transform of /i_. From Proposition 12.51 we observe that 

/—R POD 

e'^'x-^-^dx = {-lf+\Rzf / e-'^'x-^-'dx 
00 J Rz 



with the convention (—1)^+^ = e^^^'^^'^^. Therefore, R ^Im((— l)^+^a/3)//3(z) is the singular term 



contributed by yU_. Altogether, the singular part of J-"^ is i? ^ ( {Iiaaj3)fi3{z) + Im((— l)''+-'^a/3)/^(z) 
For /3 ^ N, we have 



,[/3]+l / l-(-l)/3-[/31-l\ (i?z)[^l + l-(i?^)/3 

+ 15 775 TI5-1 7^^^ r^/3 



/3---(/3-[/3] + 2) 



/3---(/3-[/3] + 2) V" /3-[/3] 
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If /3 e N, the first term is 0. Therefore, the singular terms of J^^{z) can be bounded as 

rW , m , [,3]+n , r^ B 

^ r(/3+i)'^^ ^^ ^ ^ r(/3+i)^ • 

(Step 2) The estimation of T^_ follows from the complex conjugate of ^^+, so that we only 
have to estimate the coefficient of each z^ appearing in X]«g5/3>o^'^ Ir e*^^x~^~^(ia;, excluding 
the singular terms. For each /3 > 1, first let us focus on 

w.^_/M^v (^-R^r^ .R. fry V i^R^r^^'^-' 1 

which appears in f l2.2p with additional factor (^)^. The coefficient of z"' in /i(/3) can be 
calculated as 

min{n,[/3]-2} 

r \P 



fc=0 

r\P d{f3) 



RJ ^ /3---(/3-A;)-(n-A:)! 



+ {tRY 



.Rj /3...(/3-[/5] + 2) (n-/3)-(n-[/3] + l)!' 

where 

^^^^^|l, [/3]=n + lor [/3]<n-2, 

lo, otherwise. 

We note that ,...(,_4(,_,), ^ ^ ■ ^^^^g^ ■ (^^ = ^ ■ 5(/3 + 1, fc + 1) ■ ^^, where 

S(/3 + 1, /c + 1) is the beta function f^ x'^i^l — x)^dx. The last expression can be bounded by 
just ^ ■ fcifalfcij - To estimate the coefficient of z" in the Fourier transform of /i+, we have to 

sum Coe/(Ji(/3),z") over (3 ^ S. The summation E/3g5(*^)" [ji) EfcTo^"''^^"^^ /3-(/3-fc).(n-fc)! 
can be estimated as 

\r) ~^ ^ kUn-k)\ ^ n! ' 
^es fc=o ^ ^ 

Since ^...(^_^[^]^2) " (n-/3)(n-[/3]+i)! ^ ^" (n-[^])![/3]! ' ^^^ Summation of the last term of Coe/(/i(/3), 2:") 
over B E S can be estimated as -< ^. 

"^ n! 

(Step 3) To estimate the coefficient of z'^ , we focus on 



M/'):- ,...r,l.,J --W''""'J^-)^+ E 



Z 



/3... (/?-[/?]+ 2) \ — ^ ^ ^^^^^^^^ kl{k-P + [^]-l) 

which comes from (12. 2p . For 7 G S", the coefficient of z'^ in /2(/3) is nonzero only if /3 = 7 and, 
in that case, it can be estimated as -< yt^z'^. 

Through the above steps, the desired generalized power series of J^^ has been obtained. 

(Step 4) Conversely, let us assume the conditions (1) and (2). Combining these two, we 
obtain 

l-^^WI ^ E rY^^^' ^ ^^'''' ' > 0' (2.8) 

■yes ^ '' ^ 



for constants C,B > 0. The second inequality is proved as follows. For simplicity, we suppose 
A > 1. From the condition (S2), IS* fl [n, ri + 1)| < c""'"^ for a constant c > for any n & N. 

If () < 7 < 1 then V ^'' 7'^ - T°° V ^^ --7 < V°° (cA)"+i n <- ^ /IpcAz jr 

11 u < z < i, xnen 2^^g5. r(7+i)^ - 2^n=o Z^7e.sn[n,n+i) r(7+i)^ - 2^n=o r(n+i) ^ - '^^^ • ^^ 
z > 1, we similarly get X]7g5 r(^+i) ^"^ — c^-^^^'^^^- 

Now we apply Levy's inversion formula to calculate /i, following the proof of [S]. We 
introduce an analytic function 



in C\(— oo, 0]. We also define L^{z) on (— oo, 0] to be the limit from the upper half-plane. 

For X > S and A^ e N, let f^{x) := ^^ J^ e-'''Lf,{z)dz. From the estimate ([23D, /^^(a;) 
converges to f~^{x) := ^^ S7G5 ^^ as A^ — )■ 00 locally uniformly in (S, co). By the way, chang- 
ing the path of the contour integral, one gets f^{x) = ^ Jq e~''^^J^^{z)dz — J^ e~^^L^{z)dz, 
where F^r = {A^e*^; < 9 < |}. The second integral converges to as A^ —t- 00 locally uni- 
formly. Therefore, ^ /q e~"''J-'^{z)dz — )■ 2^X^765^;^ locally uniformly in (5, 00). Summed 
with the complex conjugate, this convergence implies 



^ f e—:F,{z)dz-^-y^^ 

locally uniformly in (B,oo). Using f^(x) := j_^e~^^Lf^{z)dz, we can similarly prove 

'-:F,{z)dz^^Y.^u,{c,[^ \ 



1 
2^ 



N 
-N 



yes 

locally uniformly in (—00,— 5). Levy's inversion formula implies that fi has the absolutely 
continuous density ^ '^^es -yyo^^^ ('^7 ix) ) ^°^ l-^l > -^- ^ 

Example 2.8. (1) A Pareto distribution with the density cx~°'~^ on x > R > belongs to 
Vsa for a > 0, a ^ N, where Sa '■= {m + na : m,n G N}. Indeed, we can take Oq, = 
__c__^i(a+i)n g^j^^ Q^ _ g ^Qj. f3 ^ a in fILip . If a G N, a logarithmic term appears in the 

Fourier transform and the Pareto distribution does not belong to Vs^ (see Proposition 12.51 
and Theorem I2.6p . In addition, a Pareto distribution is infinitely divisible since its density is 
completely monotone |2U] . 

(2) As mentioned in Introduction, any a-stable distribution with < a < 1 belongs to Vs^ 
and Cauchy distributions also belong to Vn- Under the notation of (11. 7p . the density can be 
written as 

i E y"^^^^i±^^.o,f.»(i)"™'). ./O. (2.9, 

IT ^-^ mini \ \x J 1 

(m,n)GN2\{{0,0)} \ ^ ^ / 

In the case 1 < a < 2, however, an a-stable distribution does not belong to Vs^ since the double 
series f l2.9p is not convergent (the series (12.91) is still true as an asymptotic expansion [201 [2S])- 
By contrast, a-stable distributions in non-commutative probability belong to Vs^ fo^' ci-ny a G 
(0,2] (see Section E]). 

(3) Let z/ be a probability measure on [0, 00) and < a < L Then the function j^ e~'^'°^V(rfx) 
is the Fourier transform of an infinitely divisible distribution /i (see Corollary 10.6 of [22]). The 
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probability measure /x is a mixture of symmetric a-stable distributions. If, moreover, v is 
compactly supported, \x belongs to Vs^- Indeed, we have 

;-l)"m„(z/) 



^ — ^ n\ 

n=0 

and from Levy's inversion formula, the density of /i is 



This implies that yU G P^^. 

We do not treat stable distributions with index one which are not Cauchy distributions. 
Such a probability measure has a logarithm function in its density, so that it is not in the scope 

of dn]). 

2.3 Probability measures with supports bounded below 

In applications, power laws with supports bounded below form an important class. For instance, 
some stable distributions are supported on the positive real line. In this subsection, we consider 
when a probability measure /i e Vs^ has a support bounded below. 

Following the definition of Vs-, let us introduce a class of probability measures Vg whose 
supports are bounded below. 

Definition 2.9. A class Vg of probability measures consists of any probability measure \i 
satisfying the following properties: 

(1) the support of /i is bounded below; 

(2) there exist {6/3}^g5^/3>o C M and r, i? > such that 

\Ax>R{<i'^) = X] ^/3 ( ~ ) ^^' (2-10) 

|6^|<r^ (2.11) 

r < c-^R, (2.12) 

where c is the constant in the condition (S2). 

A probability distribution of V^ sometimes contains a logarithmic term in its Fourier trans- 
form, like the Pareto distribution with parameter a G N. However, a stable distribution on 
the positive real line does not have a logarithmic term in the Fourier transform. To solve this 
problem, we focus on the case S = Sa (a > 0). Then we consider when /i G Vg^ is contained 
in Vs^, which means that the Fourier transform of fi can be written as (11.41) . A crucial concept 
is a Diophantine approximation in number theory, whose importance in the context of stable 
distributions was clarified in [121 UHl US]- 

Definition 2.10. Let us consider the following condition on /3 ^ Q: for any 6 > 1, the 
inequality 

<r« (2.13) 






holds for infinitely many pairs {p,q) of Z x (N\{0}). We denote by V the set of all /3 ^ 
satisfying the above condition. 
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Remark 2.11. (1) The Lebesgue measure of V is zero from Khintchine's theorem 

(2) The set V is smaller than the set C used in [121 [IHl [13] • Indeed, C is defined as in 
Definition 12. 10[ just by replacing "for any h > 1" by "for some b > 1" . 

The following properties are basic. Let us define {(3) := min{|/3 — n| : n G Z} for /3 G R. 

Proposition 2.12. (1) Let /3 ^ Q. Then, 13 ^ V if and only if there zs A > such that 
I . / a ^ < y4" for any n G N, n > 1 . 

(2) IfxeV and z e Q\{0}, then zx,z + x eV. 

(3) IfxeV, thenx-^ G V. 

Proof. (2) can be proved from the definition of V. (3) can be proved in a way similar to [TI]]. 
Let us prove (1). If /3 ^ D, then there are 6 > 1 and go ^ f^ such that \l3q — p\ > qh~'^ for q > qo 
and p G Z. Therefore, we have (Pq) > qh~^ for q > go- This implies the existence of A > 
such that the inequality , .^/g -, < A'^ holds for any g G N\{0}. The converse is similarly 
proved. D 

The main theorem is the following. If the coefficient 6^ in fl2.10p is not zero for a /3 G N, 
then the Fourier transform of fi has a logarithmic term as in Proposition 12.51 Therefore, we 
assume that 6/3 = for any /3 G N. 

Theorem 2.13. Let fi G V§ and a ^ T). If the coefficients hp in Ii2.10\) are zero for /3 G N, 
then /i G Vs^- Conversely, if a E V, we can find a probability measure fi G Vg^ fl iVs^Y with 
b(s = for (3 eN. 

Proof. Suppose that a probability measure fi satisfies the conditions f l2.10p -f l2.12p . Then the 
coefficients {a/3}/3g5^,/3>o C C in flLip should be written as 

lmap:=bp, Rea^^ := - cot(7r/3)6/3, /3 G S'„\N. (2.14) 

The coefficients ap for /3 G N are defined to be 0. Let us suppose a ^ V. If a G Q, then 
sup{| cot{Ti(3)\ : /9 ^ N} is finite, so that a^ can be estimated as |a/3| < C^ for a constant C > 0. 
If a ^ Q U P, there is a constant C > such that |a^| < C^ since a satisfies the estimate in 
Proposition 12. 12l ( 1) . Therefore, in both cases a G Q and a ^ Q, fi can be written as 

(J-llxiyAidx) = ^ Im 

/3e5c,/3>0 

for an y4 > 0. 

Next, let us assume that a eV. For instance, let /i be defined as 

^{dx) = J2 b^x-^-^dx 

/3e5c\N 

supported on [R, oo) for some i? > and 6 > 0. The coefficients ap in f 1 1.1 1) should be written 
as fl2Aill . so that 

\map:=b^, Reo/j := - cot(7r/3)6^, [3 e Sa\fi. 

From Proposition l2. 121 (1). there is no A > such that A~l^~^\ cot{'K(3)\ is bounded for (3 G S'q,\N. 
This means that \aii\ cannot have an estimate of the form |a/3| < D^ , so that /i ^ Vs^,- □ 
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Example 2.14. (1) The a-stable distribution on [0, oo) belongs to Vg^ HP^^ for any a G (0, 1). 
If parameters in fll.6p satisfy (7, c, (3) = (0, cos(^), 1), the probability density can be calculated 
as 



ly.(-l 



1" ^ sin (7r«n)r (no; + 1) 



-l-na 



n=l 



as a convergent series for x > 

(2) (The supremum of an a-stable process [121 [ID]) Let X = {Xt)t>Q be an a-stable process 



with parameters (7, c, (3) in (11 .Op satisfying 7 = and c = w 1 + /^^ tan^(^). Let p be defined 

by p = P{Xi > 0). The process St := sup{Xs : < s < t} is called the supremum process. 
From the self-similarity of Xt, the distribution of St is the same as that of t^''°'Si. If a G (0, 1) 
and a ^ CUQ, the probabihty density of 5*1 is calculated as a convergent double series 



^P{S, <X)= X-1-" Yl ^-,n+l^""""" 

m,7i6N 



for x > 0, where bm,n is defined as 

(_l)n^+n ZL sin {^{ap + j - 1)) A sin(7ra(p + J - 1) 



T{l + ^ + n)Ti-m-an)fJ^ sin(f) fj^ sin(7raj) 

Since V C C, we conclude that the distribution of Si belongs to V§^ ^'Psc from Theorem 12. 131 
(3) (The last passage time [13]) Let Y = (Yt)t>o be a symmetric a-stable process in M'^ 
whose Fourier transform is £'[e*'^'^*] = e"*"""" for v G W^, where || ■ || is the Euclidean norm. We 
define the last passage time Ur := sup{t > : \\Yt\\ < r}. Then Uar has the same distribution 
as a'^Ur- If 1 < a < (i, the probability density of f/2 can be calculated as 



d+2m 



lp(u <t)= ^ V Jziri^zt^r 

dt ^ 2- ; ^r(^) ^,m!r(^+m + l) 

^ z ' m>0 ^ ^ ' 

which is convergent for t > 0. From Proposition 12.12( 3) and Theorem 12.131 this distribution 
belongs to P| _^ n Vs^^, iia^V. 
More examples can be found in 



2.4 Generalized moments and basic properties 

It is important to mention the uniqueness of the coefficient of z^ . 

Proposition 2.15. Let p G Vs- If J^p,{,z) = 'Ejesri^^^^^ = 'Eyes T(^)^^ ^^ f^^^ z > 0, 
then c^ = d^ for any 7 G S*. 

Proof. This can be proved inductively in terms of asymptotic behavior as 2; \ 0. D 

In the paper [9], we introduced the concept of complex moment. The class Vs enables us 
to generalize this concept. 

Definition 2.16. Let 5 be a set satisfying (SI) and (S2). The complex number c^ in Theorem 
12.61 is called a 7-complex moment of p. It is denote by m^{p). If X is a random variable with 
the distribution p, we also write rn^{p) = m^{X). 
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Remark 2.17. (1) If /i is compactly supported and 7 G N, the above 7-coniplex inoment 
coincides with the usual moment. Moreover, if /i G "Pn, ni^{^) coincides with the 7th complex 
moment [0]. Therefore, no confusion arises if we use the same symbol m^{fi) in this paper. 



From the proof of Theorem 12. 6[ we conclude the following. 

Proposition 2.18. In U.l\) . we can take a^ to be a^ = ^m^{fi). We note that z/7 G N, the 
real part of a^ is arbitrary. If •y ^N, a^ is uniquely determined by /i, so that a^ = ^m^{fi) is 
the unique choice. 

The binomial- type expansion is true for m.y{fi). 

Proposition 2.19. (1) Vs is closed under the convolution *. 
(2) For fj,,^ gVs, we have 

7,<5e5:7+<5=/3 V / ^ ; V ^ ^^ 

Proof. For ii.v e P5, let us define cp := T.^^s&S:^+5=i3 r(7+i)r(<5+i) "^t {lA'^&i.^) ■ Then, formally, 

J2fl^)(''^' = ^^('^^'^('^^ ^>0. (2.15) 

To prove that the LHS is an absolutely convergent series, we estimate C/3 as \cp\ < D^. Let 
c > be the constant in the condition (S2), i.e., \S n[n,n + 1)| < c""^^, and let A > be a 
constant such that |m-y(/^)|, |m-y(i/)| < A^. Then 

Id < V V r(/j + 1) 

'"^'-^ ^ r(7 + i)r(<5 + i) 

n=0 -y(zsn[n,n+l), ■y+5=l3 ^ ' ^ ^ ' 

<f V r(/3 + 1) 



-^ r(n + l)r([/3]-n + l) 



n=0 



= c([/3] + 2)(c+l)[^U'^. 
Therefore, f l2.15p holds as convergent series. This in addition implies that fi * u & Vs and 

3 Power laws in non- commutative probability 

We focus on free, monotone and Boolean independences which are important independences 
in non- commutative probability. In this context, the Stieltjes transform is more relevant as a 
moment-generating function than the Fourier transform. Therefore, we consider generalized 
power series expansions for Stieltjes transforms and related transforms. 
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3.1 Preliminaries 

We summarize preparatory concepts and results. Free, monotone and Boolean independences 
for random variables were introduced in [21], [16] and [21], respectively. For probability mea- 
sures n and z/ on M, the free convolution /j.Si' is defined as the distribution oi X + Y where X 
and Y are self-adjoint and free independent random variables with distributions fi and z/, respec- 
tively. The monotone convolution fi^u and the Boolean convolution /it+Ji/ are defined similarly, 
with free independence replaced by monotone and Boolean independences, respectively. 

The function -F)j(^) = (j^^ is called the reciprocal Cauchy transform of a probability mea- 
sure /i. The function F^ has a right inverse F~^ defined in F^m := {z G C_ : Im^; < 
— M, Im^; < 77|Rez|} for some t],M > 0. The Voiculescu transform 0^ is then defined as 
(p^iz) = Fij^^{z) — z. The Voiculescu transform in free probability theory plays the role of 
logJ-"^ in probability theory. The following characterizations (1), (2) and (3) were proved in 
[5], [T7] and [21j, respectively. 

Theorem 3.1. For probability measures /i and v on R, vje have the following. 

(1) (fy^mu = <Pfj. + 'Pu in an open set of the form T^j^m- 
(2)F^^, = F^oF,inC_. 

(3) F^^,{z) = F^{z) + F,{z) -zmC. 

Therefore, the transforms 0^(-2) and F^{z) — z play the roles that the logarithm of the 
Fourier transform does in probability theory. 

In free probability, stable distributions were introduced in [Sjj as analogues of stable distri- 
butions in probability theory. A probability measure fi is called ffl-stable if for any a,b > 0, 
there are c > and c? G M such that {Dafi) ffl {Dbfi) = {F>cfJ') ffl Sd- If c? is for any a, b, then 
the distribution is said to be strictly ffl-stable. In the Boolean case, l±)-stable distributions and 
strictly W-stable distributions can be defined by replacing ffl with l±) f2T\ |2]. In the monotone 
case, only strictly stable distributions have been defined [S]. A probabihty measure /x is said to 
be strictly >-stable if it is >-infinitely divisible and the corresponding >-convolution semigroup 
{fJ't}t>o with jii = n and /xq = ^o satisfies the self-similarity: for any t > 0, there is a A > 
such that /ii = -Da/^O 

The above stable distributions are characterized as follows [21 El El 1^ . 

Theorem 3.2. (1) If a probability measure fi is S-stable (resp. [^-stable), then ip^iz) (resp. 
z — F^{z)) is one of the following forms: 

(i) -7 + 6^1-'" mC_ where -i G M, 6 G C\{0},arg6G [(l-a)7r,7r] and a G (0,1). 

(ii) c — 6 log 2 in C„ where c G C_|_ U M and b > 0. 

(Hi) -7 + bz^-"" m C_ where ae^,be C\{0}, arg6 G [0, (2 - a)-n] and a G (1, 2]. 

(2) Only for strictly \>-stable laws, we change the definitions of powers: z'^ is defined to 
be e'^'°s^ in C\[0, oo) so that Im(log2;) G (— 27r,0). // a probability measure fi is strictly >- 
stable, then F^(z) = {z°' — bY^" in C_ where {a,b) satisfies one of the following conditions: 
be C\{0}, argfoG [(l-a)7r,7r] and a G (0,1];6g C\{0}, arg6 G [0, (2 - a)7r] and a G [1,2]; 
6 = 0. 

The parameter a is also called a stability index. In the case (ii), a stability index is defined 
to be one. We do not consider the case 6 > in (ii) as in the case of probability theory to 



""^In [25] J.-C. Wang defined strictly i>-stable distributions in a different way and proved that its definition is 
equivalent to the definition of [B]- 
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avoid a logarithmic term of the Stieltjes transform. A 1-strictly stable distribution is a delta 
measure or a Cauchy distribution in any case of probability theory, free, Boolean and monotone 
probability theories. 

3.2 Characterizations in terms of the Stieltjes transform and related 
transforms 

To apply the class Vs to non-commutative probability theory, we prove an analogue of Theorem 
12.61 for the Stieltjes transform and its reciprocal. 

Theorem 3.3. Let n he a probability measure and c be the constant in the condition (S2). 
Then the following are equivalent. 

(1) /i G Vs. 

(2) There are {d^)^^s C C with d^ = 1 and A > Q such that \d^\ < A^ for any 7 G S* and 
Gt,{,z) = Y.-yes ^ ^'^ {z eC-: \z\ > cA}. 

(3) There exist {h^)^^s C C (feo = 1) O'^d A > such that \h^\ < A^ and F^{z) = z^ g ^ 
in {z G C_ : \z\ > cA}. 

Proof. (1) -v^ (2): We assume that fj, G Vs- From Theorem 12.61 and Definition 12. 16[ Tp{z) has 
the generalized power series expansion J-'^{z) = ^^^zs fr'-^i'^z'^ for z > 0, and moreover we 
have |m^(/i)| < A^ for some A > 0. By the way, we can see that 



/■oo 

/ T^{z)e-y'dz = -iG^{-iy), y > 0. 
Jo 



The estimate (12. 8p is applicable to the present case, to conclude that 

M^)e~''dz = J2^^^, y>cA. 
yes y 

Since \G^{—iy) and Yli-yes "^7+1 are both analytic in {cA, cxd), we have the equality G^{z) = 

E,.5^in{^GC_:|.|>cA}. 

Let us prove the converse statement. Since the limit \miy\fiG^{x — iy) is locally uniformly 
in R\[— cA,cA], we can use the Stieltjes inversion formula, to conclude that /i is absolutely 

continuous in ]R\[— cA, cA] and ^Ji\\x\>cA{dx) = ^ X17G5 ■'■^ ("^7(/^) ij^^ ) dx. 
(2) <^ (3): taking the reciprocal of G^{z) in Theorem 13. 3[ we obtain 

y 7>0 \7>0 

It is not easy to estimate the coefficient of z^ in F^, so that we consider another proof based 
on Proposition 12.41 We can observe that this series expansion is absolutely convergent for large 

1^1 . Indeed, if \z\ is large enough, the series 1 + 'Yl,'y>o T^ + (X]7>o T^ ) + ' ' ' converges to a 
finite real number. Therefore, the order of the summands can be changed and the reordered 
series F^{z) = zJ2^(zs ^ is also absolutely convergent for large \z\. From Proposition 12. 4[ the 
coefficients b^ are bounded as \b.y\ < AC^ for an A > 0. The converse statement is similarly 
proved. D 
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The above result supports the relevance of the factor r(7 + 1) in Theorem 12.61 Indeed, the 
following is immediate, while it is not trivial a priori. 

Theorem 3.4. Let /i G Vs- If we expand G^{z) = J^jes 7^ /^^ large \z\, Im 2; < 0, then 
d^ = ?Ti^(/i). 

Remark 3.5. (1) This coincidence of m^{^) and dj supports the definition of 7-complex 
moments: one can define the same 7-complex moments both in terms of Fourier and Stieltjes 
transforms. 

(2) Theorem 12.61 and Theorem 13.31 explain the reason why G^ is defined on C_ rather than 
on the upper half-plane C+. If we stated Theorem 13.31 for G^ on C+, then the coefficient d^ 
would not be equal to m^{fi). Therefore, (0, 00) for the Fourier transform corresponds to C_ 
for the Stieltjes transform. A similar observation is in [9] in the case 5* = N. 

In free probability, the Voiculescu transform is crucial to analyze the additive free convolu- 
tion. Therefore, we characterize Vs in terms of the Voiculescu transform. For general S, the 
inverse function is difficult to treat, so that we only focus on Sai^...^ap- A key to the proof is to 
look at the powers {z"'^}^^^ as independent variables, while all of them are functions of z. 

Theorem 3.6. Let < ai < ■ ■ ■ < ap for p E N,p > 1, Sai,--,ap '■= {no + niai + ■ ■ ■ -|- UpUp : 
no, ■ ■ • , rip G N} and c be the constant in the condition (S2). Then the following conditions are 
equivalent. 

(2) There exist (e-y)^g5^^ ...^^ C C and A > such that \e^\ < A^'^'^ and(j)fj,{z) = J2-yeSc ^ 2 
in {z G C_ : \z\ > cA}. 

Proof. (1) =^ (2): let us expand G^ as in Theorem 13. 3r 2). This expansion can also be written 

as 

-1 ^ / 1 \ "" / / 1 \ "p \ "p 

(no,-,np)GNP+i 



z ^ — ' \z J \\z 



where dnQ,---,np may not be unique. In the above, formulae z'^^^ = z"'z^ and z^^ = (2;'^)" were 
used; see Remark [2.11 Let us define g{z) := G^(^) for z G C+ with small \z\. If z{l + f{z)) is 
an inverse of g{z), f satisfies 

= f{z) + J2 <,..,„,,2''°+"'°'+'"+"''°''(l + /(2))'^o+"i°i+ +"p°p+i. 

(no,---,"p)eNP+i\{0} 

To prove the existence of such an f{z), we define 

g{zo,---,z.p,w) ■.= w+ Yl <,-,np^o''^r---^;''(l + H"''+"^°^+ ^"''°''+S 

(no,-,np)eNP+i\{0} 

where (1 + w)^ is defined by the series expansion Yl'^=o (n)""^"- 9 ^^ analytic around G C^"*"^, 
and moreover, 'g{0, ■ ■ ■ , 0) = and g:^(0, ■ ■ ■ ,0) = 1. From implicit function theorem, there is 
an analytic mapping / around G C''^"*"^ such that 

gizo,--- , ZpJ{zo, ■ ■ ■ Zp)) = 0. 
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/ has an expansion of the form /(zq, ■ ■ ■ , Zp) = Y.{no,-,np)&np+^\{o} /no,-,np^o° ■ ■ ■ Zp" which is 
convergent around 0. Then we can define / to be f{z) := f{z, 2;"^, ■ ■ ■ , z'^p). 

Thus, the (right) inverse function of F^ exists as j— tttt- The Voiculescu transform (p^{z) 

is then equal to - — Trfr, which has a generahzed power series with the desired form. We note 

that the expansion of ]:.i is convergent for large \z\ with z G C_ as discussed in the proof of 

(2)^ (3) of Theorem 1331 ' 

The converse implication is similarly proved. D 

Now we have basic properties of Vs or Vs^ „ with respect to free, Boolean and monotone 
convolutions, in addition to Proposition 12.191 

Corollary 3.7. (1) For any set S with conditions (SI) and (S2), Vs is closed under l±l and >. 
(2) Vsa _....„ is closed under ffl. 

Proof. The claims for l±l and ffl are immediate from Theorem 13.3( 3) and Theorem 13. 6[ respec- 
tively. The claim for > can be proved with an argument similar to that of Theorem 13.31 For 
/x, 1/ G Vsi the function F^{Fjj{z)) has an absolutely convergent series. Therefore, we can re- 
order the summands and then F^{Fi,[z)) is of the form zY^^g^. Since this is absolutely 
convergent. Proposition 12.41 implies the existence of an A > such that j/i^l < A"' . D 

Example 3.8. (1) Stable distributions for free. Boolean and monotone independences with 
stability index a G (0, 1) U (1, 2] belong to Vs^- Cauchy distributions also belong to V-^ which 
are 1-stable distributions for any independence. 

(2) Probability distributions /i^^ defined by the Stieltjes transforms 

GUz) = r^/" ('i^ii^^iil!!^'! (, e C_) 

were investigated in [I], /i^^ is a probability measure if either of the following conditions is 
satisfied: 

(i) 1 < r < 00, < a < 1 and (1 — a)TT < argfo < vr; 

(ii) l<r<cx), l<a<2 and < argfo < (2 - a)7i. 
H'^j. G Vs^ since G'^j.{z) has a convergent series of the form ^ X^^o '^"('^' ^'^) (l) ^'^^ large 
\z\, z G C_. 



I ) 



Acknowledgements 

The author thanks Izumi Ojima and Hayato Saigo for discussions on stable distributions and 
physical applications of power laws. He also thanks Hayato Chiba for a question on a logarithmic 
singularity in the Fourier transform. Uwe Franz brought the author's attention to the relation 
between Fourier transforms and Stieltjes transforms, which worked successfully in this paper. 
Octavio Arizmendi's question on power series expansions was useful to improve this paper. 
Part of this paper was motivated by talks and discussions in 35th Conference on Stochastic 
Processes and their Applications. In particular, the author learned much from discussions 
with Juan Carlos Pardo. This work was supported by Grant-in- Aid for JSPS Research Fellows 
(21-5106). 

16 



References 

[1] O. Arizmendi and T. Hasebe, On a class of explicit Caucliy-Stieltjes transforms related 
to monotone stable and free Poisson laws, in preparation. 

[2] O. Arizmendi and T. Hasebe, Semigroups related to additive and multiplicative, free and 
Boolean convolutions, larXiv:1105.3344i ^2. 

[3] P. Bak, C Tang and K. Wiesenfeld, Self-organized criticality: an explanation of 1// 
noise, Pliys. Rev. Lett. 57 (1987), 381-384. 

[4] H. Bercovici and V. Pata, Stable laws and domains of attraction in free probability theory 
(with an appendix by Philippe Biane), Ann. of Math. (2) 149, No. 3 (1999), 1023-1060. 

[5] H. Bercovici and D. Voiculescu, Free convolution of measures with unbounded support, 
Indiana Univ. Math. J. 42, No. 3 (1993), 733-773. 

[6] U. Frisch, Turbulence: the legacy of A.N. Kolmogorov, Cambridge University Press, Cam- 
bridge, 1995. 

[7] B.V. Gnedenko and A.N. Kolmogorov, Limit Distributions for Sums of Independent Ran- 
dom Variables, Addison- Wesley Publ. Company, Inc., 1954. 

[8] T. Hasebe, Monotone convolution and monotone infinite divisibility from complex an- 
alytic viewpoint, Infin. Dimens. Anal. Quantum Probab. Relat. Top. 13, No. 1 (2010), 
111-131. 

[9] T. Hasebe, Analytic continuations of Fourier and Stieltjes transforms and generalized 
moments of probability measures, J. Theoret. Probab., to appear. 

[10] F. Hubalek and A. Kuznetsov, A convergent series representation for the density of the 
supremum of a stable process. Elect. Comm. Probab. 16 (2011), 84-95. 

[11] A.N. Kolmogorov, The local structure of turbulence in incompressible viscous fluid for 
very large Reynolds number, Dokl. Akad. Nauk. SSSR 30 (1941), 9-13 (reprinted in Proc. 
R. Soc. Lond. A 434 (1991), 9-13). 

[12] A. Kuznetsov, On extrema of stable processes, Ann. Probab. 39, No. 3 (2011), 1027-1060. 

[13] A. Kuznetsov and J.C. Pardo, Fluctuations of stable processes and exponential functionals 
of hypergeo metric Levy processes, ■arXiv:1012.0817v l. 

[14] S. Lang, Introduction to Diophantine Approximation, New expanded edition, Springer- 
Verlag, New York, 1995. 

[15] H. Maassen, Addition of freely independent random variables, J. Funct. Anal. 106 (1992), 
409-438. 

[16] N. Muraki, Monotonic independence, monotonic central limit theorem and monotonic law 
of small numbers, Infin. Dimens. Anal. Quantum Probab. Relat. Top. 4 (2001), 39-58. 

[17] N. Muraki, Monotonic convolution and monotonic Levy-Hincin formula, preprint, 2000. 

17 



[18] J.F. Navarro, CS. Frenk and S.D.M. White, A universal density profile from hierarchical 
clustering, Astrophys. J. 490 (1997), 493-508. 

[19] M.E.J. Newman, Power laws, Pareto distributions and Zipf's laws, Contemp. Phys. 46, 
No. 5 (2005), 323-351. 

[20] K. Sato, Levy Processes and Infinitely Divisible Distributions, Cambridge University 
Press, Cambridge, 1999. 

[21] R. Speicher and R. Woroudi, Boolean convolution, in Free Probability Theory, papers 
from a Workshop on Random Matrices and Operator Algebra Free Products, Toronto, 
Canada 1995, ed. D. Voiculescu, Fields Inst. Commun. 12 (Amer. Math. Soc, 1997), 
267-280. 

[22] F.W. Steutel and K. van Harn, Infinite Divisibility and Probability Distributions on the 
Real Line, Monographs and textbooks in Pure and Applied Mathematics 259, Marcel- 
Dekker, New York, 2004. 

[23] C. Tsallis, S.V.F. Levy, A.M.C. Souza and R. Maynard, Statistical-Mechanical foundation 
of the ubiquity of Levy distributions in nature, Phys. Rev. Lett. 75 (1995), 3589-3593. 

[24] D. Voiculescu, Symmetries of some reduced free product algebras. Operator algebras and 
their connections with topology and ergodic theory, Lect. Notes in Math. 1132, Springer, 
Berlin (1985), 556-588. 

[25] J.-C. Wang, Strict limit types for monotone convolution, preprint. 

[26] V.M. Zolotarev, One- dimensional Stable Distributions, Translations of Mathematical 
Monographs 65, Amer. Math. Soc. Providence, Rhode Island, 1986. 



18 



